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'^ ■ Abstract 

j!^ ■ We implement the Brink- Axel hypothesis for the excitation of the double giant dipole reso- 

^ \ nance (DGDR): The background states which couple to the one-phonon giant dipole resonance 
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are themselves capable of dipole absorption. These states (and the ones which couple to the 
two-phonon resonance) are described in terms of the Gaussian Orthogonal Ensemble of random 
matrices. We use second-order time-dependent perturbation theory and calculate analytically 
the ensemble-averaged cross section for excitation of the DGDR. Numerical calculations illumi- 
nate the mechanism and the dependence of the cross section on the various parameters of the 
theory, and are specifically performed for the reaction ^osp]^ _|_ 208p|^ ^^^ ^ projectile energy of 640 
^ ■ MeV/nucleon. We show that the contribution of the background states to the excitation of the 

DGDR is significant. We find that the width of the DGDR, the energy-integrated cross section 
and the ratio of this quantity over the energy-integrated cross section for the single giant dipole 
resonance, all agree with experiment within experimental errors. We compare our approach with 
that of Carlson et al. who have used a similar physical picture. 
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1 Introduction 

In peripheral heavy-ion colhsions at bombarding energies of several hundred MeV/nucleon and 
more, the nuclear giant resonances, in particular the isovector giant dipole resonance, are excited 
by Coulomb excitation. Recently, the excitation of higher dipole modes - the multi-phonon 
excitation - has attracted particular attention. 

In the simplest and idealized picture, the giant dipole resonance (GDR) is a one-phonon 
state described by the collective motion of all protons against all neutrons. In addition, there 
exist also higher vibrational modes (many-phonon states) that can be excited by multi-photon 
absorption. In the harmonic oscillator approximation, the excitation energy of the n-phonon 
resonance is exactly equal to n times that of the GDR, and its width is equal to n times the 
width of the GDR provided the strength distribution is approximated by a Lorentzian. 

The double giant dipole resonance (DGDR) has been observed in several nuclei: ^^^Xe [0], 
^^^Au and ^''^ Pb 0, ^ . Compared to the predictions of the harmonic picture, the measured 
cross sections for the DGDR excitation are found to be enhanced by factors ranging from 1.3 to 
2 [^ P, 1^. At the same time the experimentally determined widths of the DGDR are close to a/2 
times the width of the GDR. These discrepancies between simple-minded theoretical predictions 
and experimental results have attracted much theoretical attention. Several mechanisms have 
been studied. We mention anharmonicities of the collective Hamiltonian [Bl |9|1 and nonlinearities 



of the external field [00]. Following several earlier papers, Carlson et al. O, 12 1 have recently 



discussed the discrepancy using the Brink-Axel hypothesis. This hypothesis states that a gi- 



ant resonance is built on top of every excited nuclear state |T^. Following work by Ko [|1J], 
these authors considered the contribution to the cross section of DGDR excitation due to the 
background states which couple dynamically to the one-phonon state. It was found that this 
contribution is sizable. 

This result is interesting and calls for further study, especially since the approach of Carlson 



et al. |Tl| , uses approximations which are plausible but not based upon an expansion in terms 
of a small parameter, see Section ||. In the present paper, we apply essentially the same physical 
picture as Carlson et al. but use a formulation which allows us to derive the DGDR excitation 
cross section within perfectly controlled approximations. The resulting formula is subsequently 
evaluated numerically. Our approach makes it possible to clearly identify and calculate the 
modification of the DGDR absorption process due to the Brink-Axel hypothesis. 

The paper is organized as follows. We describe our approach in Section ^ The background 



states which couple to the collective one- and two-phonon states are complex states. These states 
are, therefore, modeled with the help of random-matrix theory [^, [l^. In Section |^, we derive 
the expression for the cross section for DGDR excitation in the framework of second-order time- 
dependent perturbation theory. The specific way in which the Brink-Axel hypothesis has been 
implemented allows for a substantial simplification of the resulting expressions (Section p.2|) . An 
analytical expression for the ensemble-averaged cross section is derived in Section^. In Section ^ 
numerical calculations are used to illuminate our results and are specifically performed for the 
reaction ^''^Pb + ^''^Pb at projectile energy 640 MeV/nucleon [0. We calculate the width of 
the DGDR, the energy-integated cross section, and the ratio of this quantity over the energy- 
integrated cross section for single dipole absorption. In Section P, we compare our approach and 
that of Carlson et al. A brief summary is given in Section |^. 

2 Hamiltonian of the Projectile 

We consider the relativistic Coulomb excitation of the DGDR of the projectile in a collision with 
a target. 

For the Hamiltonian of the projectile, we use the physical picture shown in Fig. 1. In 
the random-phase approximation (RPA), the one-phonon and the two-phonon states of the 
projectile are coherent superpositions of one-particle one-hole (Iplh) states and of two-particle 
two-hole (2p2h) states, respectively. The one-phonon state |10) is the giant dipole mode of the 
ground state |00), and the two-phonon state |20) is the giant dipole mode of the one-phonon 
state. The one-phonon state |10) is dynamically coupled by the nuclear Hamiltonian to more 
complex particle-hole configurations \Qk) with k = 1, . . . ,K and K ^ 1. This coupling causes 
the giant dipole resonance to acquire a spreading width T{q. The two-phonon state |20) is 
likewise coupled to such configurations. In order to accommodate the Brink-Axel hypothesis, 
we group these configurations into two classes. States in the first class are denoted by \lk) with 
k = 1, . . . ,K. Each such state represents the dipole mode of the lower state \0k). States in the 
second class are denoted by |0a) with a = 1, . . . ,M and M ^ 1. These latter states are not 
coupled by the dipole operator to lower-lying configurations. The states |20), |1A;) and |0q:) are 
all dynamically coupled to each other. 
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Figure 1: Schematic illustration of the DGDR excitation and the couplings between the phonon 
states and many-particle many-hole configurations. The level spacings of the configurations are 
exaggerated in the figure. 4 



According to this picture, the Hamiltonian matrix H has the form 
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The matrix H consists of three diagonal blocks. The first block has dimension one and contains 
the energy Eq of the ground state |00). The second block has dimension 1 + K. It contains the 
unperturbed energy Ei of the one-phonon state |10), the elements W^i with k,l = 1, . . . ,K oi 
the Hamiltonian matrix /J'^°) governing the states |0A;), and the real coupling matrix elements Vik 
connecting the one-phonon state with the states |0A;). The third block has dimension 1 + K + M. 
It contains the unperturbed energy E2 of the two-phonon state |20), the Hamiltonian matrices 
H\,i with k,l = 1, . . . ,K and TiJ^j with a, (3 = 1, . . . , M governing the states |1A;) and |0q;), 
respectively, and the real coupling matrix elements V2A;, V2a and Wka connecting the three sets 
of states. In this simple picture, we pay no attention to spin and isospin. Inclusion of these 
quantum numbers is not difficult but requires a straightforward extension of our formalism. If, 
for instance, the ground state has spin and parity 0'^, then all states in block two have spin 
and parity 1~, while the states in block three decay into three groups having spin/parity values 
0+, 1+ and 2"*", respectively. 

We turn to the statistical assumptions. We describe the K states |0A;) with K -^ 00 in terms 
of the Gaussian Orthogonal Ensemble of random matrices (GOE), so that H^^' represents this 
ensemble. For excitation energies of the giant dipole resonance (which lie typically between 8 
and 15 MeV), this description seems eminently reasonable, except perhaps for the lightest nuclei. 
The spreading width V\q = 27cv'^/d of the giant dipole resonance is expressed in terms of the 
mean level spacing d of the states |0A;) and of the mean square coupling matrix element Vif. = v"^ . 
It was mentioned above that for each value of /c = 1, . . . , -fr, the state \lk) is assumed to be the 
giant dipole mode of the state |0A;). This implies that the matrix H^^^ is identical to H^'^'> except 
that the center of the semicircle is shifted by E2 — Ei, so that 



Hi]^ = {E2-Ei)6M + H, 



(0) 
kl 



(2) 



We also assume that Vik = V2k- This assumption is not strictly implied by the Brink- Axel 
hypothesis because the states |10) and |20) (or |0A;) and |1A;), respectively) are not identical. 



Rather, we have |20) oc -D|10) where D is the dipole operator, and correspondingly for the 
states labelled k. In practice, the assumption Vik = V2k implies that the spreading widths for 
the couplings |10) <^ |0A;) and |20) -v^ \lk) are equal. This assumption seems plausible. We 
return to this point in Section ^. We note that with these assumptions, the second block of the 
matrix (0) differs from the corresponding part of the third block only by {E2 — Ei) times the 
unit matrix. The matrix Ti'^^^ is also assumed to represent a GOE, with M taken to infinity. We 
assume that 7i^^^ and H^^^ are uncorrelated. 

The states in the second and third block can decay by particle emission. We take account of 
this fact by introducing the decay widths FJq and r2o of the one-phonon and the two-phonon 
states, respectively, and the decay widths FJ and T2 of the states |Ofc) and \ik), respectively. 
Within the statistical model, it is obviously adequate to assume that FJ and Fg are independent 
of the running index k = 1, . . . ,K. In Section |3.2| it is shown that it is not necessary to assign 
a decay width to the states |0a). The width matrix S accordingly has the form 
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The effective Hamiltonian Hc^ is then given by 
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All four decay widths can be calculated with the help of the optical model. We return to this 
point in Section |. 

3 Second— order Time— dependent Perturbation Theory 



We suppress the intrinsic structure of the target and replace it by a point source of the electro- 
magnetic field. To describe the DGDR excitation of the projectile, we use the standard approach 



to Coulomb excitation [|1^: The relative motion of projectile and target is described classically, 
and the intrinsic excitation of the projectile is treated quantum-mechanically. 



3.1 General Approach 

The time-dependent Schrodinger equation 

ih^^^ = (H^s + H^mm) (5) 

for the internal wave function \ip{t)) of the projectile contains the interaction Hi{t) with the 
time-dependent electromagnetic field caused by the relative motion. We consider only dipole 
excitation and write Hi{t) = Dh{t) as the product of the dipole operator D and of a time- 
dependent function h{t). We have h{t) -^ for t -^ ±oo. A suitable form for h{t) is given 
in Section |^. The dipole operator D induces transitions between the following pairs of states: 
|00) -* |10), |10) -^ |20), and |0A;) -* \lk) for all k = l,...,K. 
We use the interaction representation 

\<f){t))=exp{iH,st/n)\tlj{t)) (6) 

so that 

m^m = mn)m) (7) 

with 

Hi{t) = exp{iH,s t/h)Hi{t) exp(-zi^eff t/k) . (8) 

Excitation of the DGDR is a two-step process. We use second-order time-dependent per- 
turbation theory. With |0(— oo)) = |00) and Hi{t) — ^ for t — ^ ±oo, we have 

1 roo rti 

l0(+oo)) = i-^f / dhH,{t^) / dhHr{t2)m . (9) 

To calculate the intensity I2 for double dipole excitation with energy transfer £", we first consider 
the case where r2o and Tl both vanish. Let \n) represent an eigenstate with energy En of the 
submatrix Hcs,33 forming the third block of the matrix Hes- Then I2 is given by 

h{E' + Eo) = Y.\{n\<P{+oo))\'6{E' - {E^ - E,)) . (10) 

n 

For I2 we have used the argument E' + Eq rather than E' because it is convenient to introduce 
E = E' + Eq. Moreover, we use the identity 

T.\n){n\S{E-En)=^-^lm \ , (11) 

n Tl Hi Xleflf 33 



the definition of Hi{t) given above, and tlie fact that Hcs\00) = Eo\00) to rewrite hiE) in the 
form 

\ P OO /'Li P OO P~t 

h{E) = -J dtj dt2 dt[ / ' dt'^{00\H*{Q exp{-t{H:^ - Eo)(4 - t[)/h} 

71 ll J —OO J —OO J —OO J —OO 

xH*it[) exp{-z(iJ:ff - Eo)t[/h}lmi-—-^ ) exp{t{H,s - Eo)h/h} 

xHiih) exp{t{H,s - Eo){t2 - h)lh}H^{h)m . (12) 

The form of Eq. (|l^) remains unchanged when we take account of the finite values of the decay 
widths Fjo and Fg. Substitution of Hi{t) = Dh{t) yields 

1 /"OO /-O roo rO 

h{E) = -J / rfti / dT2 / dt\ / dr'^ h{h)h{T2 + h)h*{t[)h*{T^ + t[) 

Vr Ti J—OO J —OO J—OO J —OO 

x(00|Dexp{-2(if:g,22 - E,)T'Jh}Dexp{i{E - Eo){h - t[)/n} 

xlm(-— ^-- )Dexp{z(/Jeff,22 - Eo)T2/h}D\00), (13) 

^^ — -ncff,33 

with r2 = t2—ti and Tg = t'2—t[. We have replaced exp{i(ifeff,33 — -E'o)(ii —i'l)/^} by exp{2(£' — 
Eo){ti — t[)/h}. The remaining two exponentials in Eq. ([TsD contain the second block Hes,22 of 
Hes and can be expressed in terms of Green functions. The integrand of 

roo I 

/ de- --exp(-zet/^) (14) 

has poles in the lower half of the e plane. We may close the contour of integration in the upper 
(lower) half plane if t <0 (t >0, respectively). Hence, 

/•OO I 

/ de— -—exp{iet/h) = -2i7ie{-t)exp{iH^st/h) , (15) 

J—oo e — ^eS 
poo I 

/ de— -— exp(-ie t/h) = 2z7r0(-t) exp{-iHe{it/h) . (16) 

J—oo e — J^cS 

We observe that both the r2-integration and the r2-integration are restricted to values <0. We 
use the formulas (^5|) and (pJ| ) to find 

1 roo rO /-OO ^0 

hiE) = - (—4) / dh / dT2 / dt[ / dr^ h{h) h{r2 + h) h*{t[) h*{r^ + t[) 

TT Tl J—oo J—oo J—oo J—oo 

/OO roo 

de / de' exp{i(e - ^0)^-2/^} exp{-2(e' - Eoy^/K] 
-00 J—oo 

xM{E,e,e'). (17) 

Here, M is a stochastic function of H^s defined by 

M{E, 6, e') = mo ^ Dlm{—-^ )^;^4? ^10°) • (l^) 

l^ ) ~ ^eS,22 ^ ~ J^eS,33 e^ " -neff,22 



To simplify M, we observe that D does not connect to the rows and columns labelled a and (3 
in the third block of the matrix in Eq. (|l]). Therefore, we restrict attention in the third block to 
the subspace of states obtained by excluding these rows and columns. The projection of i/cff,33 
onto this subspace (which we denote by s) is written as Hes,s, and we have 

1 , 1 



where W = {V2a, W^'a)- ^^ insert this result back into Eq. ([T8| ) and find 



(19) 



M{E, e, e') = (00|Z^-^^-^Mm(-^— ^^, , (w'Y ^^ r R ^'°°^ " ^^^^ 

e -^efr,22 E - JJeff.s " W e+-H(") ^^ > e " ^efr,22 

Eqs. ( P^ and ( PP] ) are the central results of this Section. 

3.2 Simplification of M 

The expression for M can be simplified further because H^^^ and T-i^^^ are uncorrelated. Since 
Ti.^^^ appears in the intensity h^E) only via the expression (pil|), we can perform the average 
over Ti.^^^ right away. We have 

Here r2o is the spreading width of the two-phonon state due to its coupling to the states |0a), 
and Eg is the common spreading width of the states |1A;) due to their coupling to the states 
|0a). We observe that the result of the averaging procedure is independent of whether or not we 
include a decay width for the states |0a), see the end of Section 0. We use the result of Eq. (pi|) 
in Eq. (pO|). To justify this step, we observe that a series expansion of the right-hand side of 
Eq. (^) in powers of the Green function of T^^^-*, followed by averaging over 'H^^\ would amount 
to taking the average individually over each Green function appearing in the expansion. This is 
because the energy argument in all these functions has a positive imaginary part. Resummation 
of the result amounts to using the result of Eq. (|2T|) in Eq. (20). This yields for M 



W' ^ ^ ^,,J Wy = -Ei = -m) ^-'" " , I . (21) 




M{E, e, e') = {OOlD ^ ^^ Dlm{ J )D i D|00) . (22) 

^ ^ -"off,22 E - Hcs,s + ^s e - -ncff,22 

A further simplification of M is possible because the real part of -f/^eflF,s differs from the real part 
of Hes,22 only by {E2 — Ei) times the unit matrix, see Eq. (^. To use this fact, it is necessary 

9 



to write the matrices appearing in the three Green functions in Eq. ( ^21) more exphcitly. We 
denote the real part of Hcs,22 by H22 and have 

H^ = I ^: !;' I . (23) 



Moreover, we define 






rio 



Si = (V2) I ': ^ „, I , (24) 



hA 



and, with 120 = rl + r^, r2 = r| + vi 




(25) 



Then, M takes the form 



e — ri22 — ^1 h, — rj2 + rLii — n.22 + 1^2 e — -H22 + ^1 

A further simphfication of this expression is possible with the help of the Brink- Axel hypothesis 
and of algebraic identities. In the spirit of the Brink- Axel hypothesis, we assume that the dipole 
transitions |00) —>■ |10) and |0A;) —>■ \lk) all have identical dipole matrix elements. We denote all 
these identical matrix elements also by D. We also take account of the fact that in the harmonic 
oscillator picture, the dipole matrix element for the transition |10) -^ |20) is given by v2-D. 
Using these facts, we can express M completely in terms of the Green functions 

G±(B) = (10|— — i— — 110>, (27) 

To this end, we introduce the (1 + i^)-dimensional unit matrix X, define 6a,6ao,6T,6To,6p 
and 6pQ by 

S2-S1 = 5ct J + 5cto |10)(10|, (28) 

E2 + Si = 5r J + (5ro |10)(10|, (29) 

2Si=5pJ + 5po |10)(10| , (30) 

and make use of the identities 

1 1 1,1 1 . 



E-H22 + ^2 e-H22 + ^i E-e + 5a^e-H22 + T.i E- H, 



22 



'"° ' mm J , . , (31) 



E-t + 5a E-H22 + T.2' " ' e - if22 + Si 

10 



I I =^ ( I I ) 

E- H22 + E2 6' - i?,2 - El E-e' + St\' -H22-T.1 E- H22 + E2 

^Tn 1 , , , , 1 , ^ 

-p TT^ 10 10^— — , 32 

1 1 1.1 1 . 



e' - H22 - Si e - if22 + Si e - e' + (5p e' - H22 - Si e - i/22 + Si 

^'' ' |10)(10| ' ,^ • (33) 



t- t' + 5p t' - H22 - Si e - if22 + Si 

We define E" = E — E2 + Ei and obtain as a result 



G+(6) ,1 1 , 



e-e' + 6p^E" -e + 6a E'-e-dr' 

^^^^' V '^' '1 i~^2 {E ) — —f, — — —f, , r- G2 {E )) 

E — e — OT E — e — da E — e + da E — e + or 

, ^>o , /7T ,, G2iE")Gtie) , 5to ^^ G^ie)Gt{E") 

~^T^' A~ + V 2 - 1) —r, — J (^7^ r—^ + 1 - V2) —-j-, — — 

E — e — OT E — e — oa E — e + or E — e + da 

, ^^0 /7^ ^, G2{E")G-,{e) , 5a, ^, Gt{e)Gt{E") 

^E"-e'-5a^^ ^ E"-e'-5r ^E"-e + 6a^ ^ ' E"-e + 6T 

+G^{e')Gt{e) {G,iE")[-- ^"°^"° 



{E" -e-dr) {E" - e' - da) 
drodao 



E — e — dr E — e — da 

+Gt{E")[- 



{E" -e + da) {E" - e' + dr) 
+(^^-l)(^-^^-l-^^) + 2v/2-31}). (34) 

Terms carrying the factor (v2 — 1) arise because not all dipole matrix elements are equal to D. 
We note that the stochastic matrix iJ22 appears in Eq. (p^ ) only in the Green functions Gi and 



4 Ensemble Average and Energy Integration 

The ensemble average over H^^^ affects only the Green functions Gi and G2 in Eq. (^41). The 
terms in the first three lines of Eq. (|34D contain a single Green function (a "one-point function") 

11 



only. The next four lines contain products of two Green functions of the type G^G^ (a "two- 
point function"), and the last four lines contain two three-point functions. Calculating the 
ensemble average of the one-point functions is straightforward, see Eq. (PT]), and yields 



1 



Gf (e) = , (35) 

e — £/i ± ^i 10 



G^{E) = . (36) 

_t/ — _C/1 ± 2! 20 

Here Tio = TJo + T{^ and Tao = f^ + A^. 

We turn to the ensemble average of the two-point functions. These are written as the 
sum of two terms, the product of the averages of the two Green functions and the average 



over the product of the fluctuating parts 5G. For example, Gi{t)G2{E") = Gi{e) G^(i?") + 
6Gf{e)6G2{E"). The last term can be calculated with the help of the supersymmetry technique 
[|l^ , |I9|1 . We do not give the details here because application of the technique shows that these 
terms are negligible. This is because the decay widths FJq, FJ, r2o, r2 turn out to be orders of 
magnitude larger than the mean level spacing d of the states \0k), see Section ^. This fact 
leads to an exponential suppression of the fluctuating part. As a result we have, to a very good 
approximation. 



Gt{e)G2{E") = Gt{e) G^^E") . (37) 

We have already remarked that the average of the product of two retarded Green functions (or 
of two advanced Green functions) is trivially equal to the product of the averages. Thus, we 
find that the averages over all two-point functions in Eq. (|3^ are equal to the product of the 
averages of the two Green functions. We turn to the average over the three-point functions. 
Here, a similar argument is expected to apply. Unfortunately, the supersymmetry technique is 
not capable of giving an analytical result in this case, and we have used numerical simulations 
to estimate the fluctuating part. For the nucleus ^°'^Pb, for instance, we have found that the 
contribution of the fluctuating part is less than one part in 10^. Hence, we have 



/^+/^-/^± r^ n+n-^ n^ n+ n^ n^ ('i<A\ 

Ut]^ kj^ Ur2 — yJ\ yjr^ (J2 — >-ri ^\ 1-^2 ; V"^"J 

and the ensemble average of M{E, e, e') and I2{E) can be computed in terms of the ensemble- 
averaged one-point functions. 
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4.1 Total Intensity 

Using contour integration, we calculate the integrals over the variables e and e' . This calculation 
is lengthy but straightforward. Using Eq. (^If), we find 

n /•+00 /-O /■+00 ^0 

hiE) = 271 i{—)^ / dti / dT2 / dt[ / dr^ 

hih)h{h + r2)h*{t[)h*{t[ + r^) exp[z(E - Eo)ih - t'l)/^] 

X (/oe(r2 - r^) exp[-2(Ei - Eq + '-TioKM exp[t{E^ ~ Eq + '-T^o - i^p)r2/Ti] 

xifnE)-f2-{E)) 
+/oe(r2 - r^) expH(i^" - Eq + ^Ar)r^/;i] exp[2(E" - Eq + iAa)r2/h]{f+{E) - f+iE)) 
-/oe(r^ - r^) exp[-2(i5;" - i^o - ^Aa)r^/;i] exp[2(i^" - Eo - i/\T)r2/n]U^{E) - f^{E)) 
-fo&{r^ - ra) exp[-2(Ei - Eq - ^Fio + iAp)T^/h] exp[i{E^ - Eq - '^T,o)T2/h] 

x{fi{E)-fnE)) 



with 



+ (V2 - 1 + zAaoUiEm-G^iE" 



+ (1 + zAroG2-(^"))/r(^)) exp[-z(Ei - ^o + 2^10)^^/^] exp[^(E - Eq - 2Ar)r2//i] 
+ (1 - V2 + 2Aao/2+(i?))(-G+Cn 

2 . 



(1 - iAToGt{E"))f+{E)) exp[-t{E -E^ + Ar)T'Jh] exp[z(^i - E^ - -Tio)r2/n] 
-[(1 - v^ + zAao f+{E)){l-iAr,GtW))ft{E) 



iV2-l+iAa, f2iE))il + tAroG2iE"))f{iE)) 



1 



+Apo:p ^iUiE) - MiE)) + (V2 - 1){G2{E") - G+(E"))] 

Tio - Ap 



X expH(Ei - Eo + ^rio)r^/^] exp[z{Ei - Eo - ^rio)r2/^]) , (39) 



fHE) ^ 



E-E2T frio ± iAt 



fHE) ^ 



E-E2± frio ± lAa ' 



Eq. (|39D constitutes the main result of the theoretical part of this paper. Under perfectly 
controlled approximations, we have derived an expression for the intensity which embodies the 
Brink-Axel hypothesis and which describes the formation of the DGDR as a transport process. 
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4.2 Intraband Intensity 



We now calculate the intensity If^^^i^E) which would result if only the ground state and the 
one-phonon state could absorb dipole radiation. In other words, we suppress dipole absorption 
by the states labelled |0A;), although we do keep the dynamical coupling of all the states as 
described by the matrix ifcff defined in Eq. (^. We do so in order to distinguish the dipole 
excitation taken into account by the usual approach to the problem, from the transport process 
described by Eq. (|39D . This intensity is defined as 

2 /•+00 (-0 /•+00 (-0 



V /■+00 /-u /■+00 /-u 

/i-t-(E) = _ / dtj drj dt[ dr^h{t^)h{h + r,)h*it[)h*{t[ + T^] 

TCTl J —oo J —oo J —CO J —oo 

/OO /"OO 

de / de'exp{i{e - Eo)T2/h} 
-oo J —oo 



X exp{-z(e' - EoK/h}{00\D\lO){lO\^— — |10)(10|D|20) 

e — ii22 ~ ^1 

x(20|Im(— -— — ^-— -— — )|20) 

£/ — i!/2 + -C/l — -"22 + ^2 

x(20|D|10)(10| \ |10)(10|D|00) . (41) 

The ensemble average is given by 

D p+oo pO p+oo 



/ ) r+oo pU n+oo pU 

/intra(^) = 277 X 2 X ( — )^ / dti / ^Ta / dt[ / dr^ h{ti)h{ti + T2)h*{t[)h*{t[ + r^) 

ri J— oo J — oo J —CO J —OO 

X exp[z(E - Eo)ih - t[)/h]jT^ J'V iF2 exp{z(Ei - Eo - t/2T,,)T2/h} 

[E - ii2J + if 20 

X exp{-z(Ei - Eo + z/2rio)r^//i} . (42) 

The factor 2 arises because the dipole transition |10) — ^ |20) is twice as strong as the transition 
|00) -^ |10). We note that when Fio and r2o go to zero, J™*''^(-E) is identical to the harmonic 
limit denoted by /2^''(E). 

For the discussion of our results and the comparison with experimental data, it is useful to 



also give the average intensity Ii{E) for the excitation of the GDR. It is calculated similarly and 
given by 



1)2 ^+oo p+oo p 

2^J-oo "^^'J-oo ^^1 Hti)h*{t[) exp[z{E - Eo){h - t[)/n] ^^ _ ^^|° ^ ,^,^ 



4.3 Cross Sections 

The cross section for the DGDR excitation is given by 

POO 

a2{E) = 27i hdhhiE). (44) 

*^t)niin 
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Here femin is the minimal impact parameter which is introduced in order to account for the strong 
absorption that occurs as soon as the colhding nuclei come within reach of their nuclear forces. 
The cross section cri{E) for single dipole excitation is defined analogously. We also define 

a2(E)'°*''^ = 271 bdb If^'\E) (45) 

•^Oniin 
/•CXD 

^^'■(^) = 27r / hdh I^^'iE) . (46) 



and 



We define two enhancement factors. The first one compares our result with the harmonic limit 

and is defined as 

^ JZo dEa,{E) 

" IZo dEar{E) • ^^'^ 

The second factor compares our result with the cross section calculated without the Brink-Axel 

hypothesis and is defined as 

_ J-^ iEa,(E) 

r„ dEafiE) ^ ' 

It is also of interest to compare a™^'^^{E) with the harmonic approximation. This is accomplished 
by the third factor 

I-oo dEa^iE) ■ ^^'^ 

The contribution of the extraband excitation is measured by the following two ratios, 

pextra _ IZo dEME) - a^r%E)) 
' ~ iZodEa^iE) ' ^'"^ 

pe.tra ^ I-oo dEja^jE) - at'-^E)) 

iZo dEa^^E) ■ ^' > 

We note that Rl^^"^^ and Rf^^"^^ also account for the contribution of interference terms between 
the intraband and extraband excitations. 



5 Numerical Results 

For the calculation of the decay widths r|o,r|,r2o and Vl, we have used the computer code 
developed by E. Sheldon and V. C. Rogers [^. It contains a global optical-model potential 
to compute the transmission coefficients of nucleons. These in turn were used to determine the 



mean absorption cross section using the exciton model |21]: The average decay rate of a state 
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with n excitons into a channel where the nucleon has (asymptotic) kinetic energy e^ is given by 
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2mek , pn^i{p-l,h,U -Eb-ek) ,_^, 

where m is the reduced mass, U the excitation energy of the residual nucleus, Ei, the separation 
energy of particle b and a{ek) the mean absorption cross section. The function pn{Pi h, U) is the 
ra-exciton state density with excitation energy U, 

^"'^•"•^) = p./,.(p+\-i). <^«)'"'"- (^^) 

Here g is single-particle level density, p and h are the numbers of particles and holes, and 
n = p + h. The decay width is given by 

r^ = nJ2 f^'"'' Wr.{ek)dek . (54) 

The decay widths rjg and FJ are evaluated at energy Ei, and F^o and F2 at energy £^2- 

In the long-wavelength approximation [^ where the impact parameter b is large compared 
to the nuclear radius r, the time-dependent function h(t) is given by 



Here r = 7 t> t/b with 7 = 1/J{1 — v'^/c^), v is the relative velocity, and uj = {Ef — E.j)/h. The 
validity of the long-wavelength approximation was checked ^3[ for giant resonance excitation 
in the process ^°^Pb + ^'^'^Pb at 640 MeV/A and was found to be accurate to within a few 
percent. 

We consider the DGDR excitation of a ^°^Pb projectile, incident on a ^''^Pb target. This 
reaction has been studied at 640 MeV/A at the GSI/SIS, Darmstadt |^. We apply the formalism 
developed in the preceding Sections to calculate the cross section, enhancement factors, and the 
DGDR width. Using the method described above and with p = h the values p = I for FJq, p = 2 
for FJ, p = 2 for F20, and p = 3 for Fj, we find for the decay widths the values FJq = 0.11 MeV, 
F|o = 0.026 MeV, F| = 0.30 MeV and F^ = 0.16 MeV. We set Ei-Eo and E2-E0 equal to their 
experimental values 13.5 MeV and 27.0 MeV, respectively. The parameter u appearing in the 
parametrization (|55D of h{t) was accordingly chosen as hu = 13.5 MeV. Various parametrizations 
of bmin with regard to the nuclear system have been proposed. Most widely used are those of 



Ref. 24 



1.34(v4i/' + Al^^ - 0.75(^7'/' + A^'/')) (fm) (56) 
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and of Ref . |p5| , 



brmn = 1.1(^1/' + A^ + f/g ^^,3 - 1.9) (fill) . (57) 

AY^ + Ar 

For the system ^^'^Pb + ^^^Pb, Eqs. (|6]) and {^) yield bmin = 15.7 fm and bmin = 14.4 fm, 
respectively. We have used the mean of these estimated values, bmin = 15.05 fm. 

In order to calculate the absolute values of the cross sections we evalute the dipole matrix 
elements by means of the sum rule for dipole transition 

2_^ujfi\Vj. I - -Z-— — - bD ■ (oSj 



r Sttttin a 



Jm) 



Here Dj^' = J rYijn{f)p{r)dv. In practice, we use 

j:ujMD'\ = ^Sn, (59) 

/ ^ 

with D = J xp{r)dv = J zp{r)dv. Unless otherwise stated, we always use 100% of the sum rule. 

Particular attention must be paid to the spreading widths as these determine largely the 
outcome of the calculation. Here, we repeat a statement already made in the Introduction: This 
paper is primarily devoted to the implementation of the Brink- Axel hypothesis, to the derivation 
of an expression for the average cross section based on this hypothesis and on controlled approx- 
imations, and to the investigation of the consequences of this hypothesis for the excitation cross 
section of the DGDR. In this sense, we consider the various spreading widths as parameters. 
We are not primarily concerned with assigning physically realistic values to the F^^'s. Still, it is 
worthwhile to comment on some problems that arise if one tries to do so. 

The parameters are the spreading width F2 of the states |l/c), the spreading width F|q of the 
one-phonon state |10) (by assumption this width also describes the mixing of the two-phonon 
state |20) with the states |1A;)), and the spreading width F20 which describes the mixing of 
the two-phonon state |20) with the states |0a). The total spreading width of the two-phonon 
state is accordingly given by Ti^^^i 20 = F|q + Fao- We know experimentally the spreading width 
Fio = 4.0 MeV of the one-phonon state. Theoretically Fio is given by Fio = F|q + F|q. 
Typically, three effects contribute to the observed spreading width of the GDR: The coupling to 
the complex modes of excitation contained in our Hamiltonian of Eq. (|l|), the spread in energy of 
the single-particle modes which contribute to the GDR (Landau damping), and the coupling to 



low-lying vibrational modes ||26|. This makes it difficult to estimate the spreading width T^^^^^ 
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of the two-phonon state. Neglecting the contribution of Landau damping and of the vibrational 
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modes to FJq altogether and using the exciton model in the form of Ref. [^ (we identify the 



n-phonon state as a ra-particle ra-hole configuration), we find that the spreading width Ti^^^i „q 
of the ra-phonon states is proportional to 2n Im [Vopt{En/2n)]. Here Vopt{E) is the depth of the 
optical-model potential at energy E, and En is the excitation energy. This relation derives from 
the fact that each exciton (particle or hole) carries the average energy {En/2n). We observe that 
for n = 1 and n = 2, the energy arguments of Vopt{En/2n) coincide, yielding Ti^^^i 20 = 2rJQ. 
This result is in keeping with the predictions of the harmonic picture for the n-phonon states 
and gives r2o = TJo = 4.0 MeV. Comparison with experimental data will show that this is an 
overestimate which we ascribe to the fact that Landau damping and coupling to the vibrational 
modes contribute differently to FJo and to Fgo- It turns out that r2o ~ 0.5 FJq yields results 
which are in good agreement with the data. Similarly, using the measured value for FJq and the 
exciton model to estimate Fg yields an unreasonably large value. We often use F2 = 0.5 MeV 
and remark only that our results are insensitive to this choice, see Fig. 5. 

Prior to our giving detailed results for the reaction ^°^Pb on ^°^Pb at 640 MeV/A (the 
case studied experimentally), we discuss the dependence of our results on some of the relevant 
parameters. The following calculations were all done for the reaction ^osp]^ qj^ ^''^Pb as described 
above except that the energy, the decay widths and the spreading widths were varied. 

We investigate the dependence of the cross sections for the GDR and DGDR excitation on 
the spreading width FJq and on the projectile energy Ep/A. It is instructive to begin with the 
GDR excitation. Fig. 2 shows the energy-integrated cross section o"i for the GDR excitation as 
a function of the spreading with FJo at various projectile energies. The dependence of ai on FJo 
changes significantly with projectile energy. For the lower projectile energies, the cross section 
first increases strongly then more slowly and finally decreases as FJo increases. For large Ep/A 
(10000 MeV) the cross section decreases monotonically with increasing Pfg. This dependence 
can be understood qualitatively with the help of the Weizsaecker- Williams (WW) method of 
equivalent photons [0, ^]. In this method the transition probability for the GDR excitation at 
impact parameter b is given by 

PcDRib) = J N{b, E^)aphot(i5,)-^ . (60) 

Here N{b, E^) denotes the number of equivalent photons of frequency E^/h at impact parameter 
6, a function of projectile energy, while ap\iot{,E^) is the photoabsorption cross section which 
is well described by a Lorentzian. For small values of the spreading width F|q, only the pho- 
tons with frequency close to the central frequency {Ei — Eq)/Ti substantially contribute to the 



transition probability. As FJq increases, an ever wider band of photons contributes to and en- 
hances -Pgdr(^)- On the other hand, the widening of the normahzed Lorentzian causes the 
absorption strength of the photons with frequencies around the central frequency to decrease, 
which suppresses -Pgdr(&)- These effects compete against each other. The result depends on the 
distribution of available photons over the relevant frequency range. For small values of Ep/A, 
N{b, E^) decreases rapidly with increasing E^ [^]. In this case the gain of Pgdr(&) exceeds the 
loss unless FJq becomes large. As a result, the GDR cross section (an integral of -Pgdr(&) over 
impact parameter) increases with F Jq up to a certain value of the spreading width beyond which 
the loss surpasses the gain. At very high projectile energy the curve of N{h, E^) versus E^ be- 



comes flat [^. The loss always exceeds the gain, and the cross section decreases monotonically 



with F{q. Since iV(6, E^) increases with Ep/A, so do both Pgdr(&) and cross section. Fig. 2 
shows that the enhancement of the cross section due to the spreading width is most pronounced 
when Ep/A is around several hundred MeV. We refer to this enhancement as to the damping 
enhancement. Taking for projectile energy and spreading width the experimental values Ep/A 
= 640 MeV and FJq = 4.0 MeV, we find for the ratio of the energy-integrated cross section with 
damping to that without damping the value 1.20. 

Fig. 3 shows for several projectile energies the dependence of the energy-integrated cross 
section for the DGDR excitation on the spreading width FJq. We use the five ratios defined 
in Section 4.3 to measure the cross sections. We keep the decay widths and the ratio FI0/F20 
fixed. We focus first on the intraband cross section measured by i?™*'''^. The dependence of the 
intraband cross section on spreading width and projectile energy is similar to that of ai (Fig. 2): 
At small projectile energies, R^^^^^ first increases then decreases with increasing Fjg, while it 
decreases monotonically with F|q at high energies. This is because the intraband excitation is 
governed by two transition probabilities: one for the transition |00) — > |10), the other for the 
transition |10) ^ |20). Each of the two probabilities can be described by Eq. (|60l). Thus, the 
intraband cross section depends on spreading width and projectile energy in a similar way as cti. 
Some of the probability flux feeding the one-phonon state |10) flows into the states \0k), and 
- in the intraband approximation - is not available for the second dipole excitation. This flow 
increases with T{q. Therefore, the enhancement of i?™*''^ only survives for rather low projectile 
energies, and the range in which _R"^t''^ increases with the spreading width is significantly reduced 
compared to that of cxi. The extraband excitation measured by i?°^*^^ (compared with the 
harmonic limit) and Rf^^"^^ (compared with the intraband cross section) increases monotonically 
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Figure 2: Energy-integrated cross section for the GDR excitation as a function of spreading 

I T 

width rjo at different projectile energies Ep/A, with Fjo = 0.11 MeV. 



20 



with rjo- For small projectile energy [Ep/A = 200 MeV) and small values of FIq both ratios 
have small negative values. However, in most cases the extraband excitation gives a positive 
contribution to the total cross section. 

The total cross section is measured by -Ri (compared with the harmonic limit) and by R2 
(compared with the intraband cross section). The ratio R2 increases monotonically with r|o and 
is larger than unity except for small projectile energy and small FJq. At the smaller projectile 
energies, Ri first increases and later decreases with Fjg. The turning point shifts towards larger 
values of F|q as Ep/A decreases. With increasing Ep/A the range where i?i is larger than unity 
disappears. An enhancement of the total cross section compared with the harmonic limit can, 
therefore, be expected only for low projectile energy. Moreover, this enhancement is not mainly 
due to the extraband excitation [Rf^^"^^ and R'^^^^ being small), but is mainly due to the damping 
enhancement discussed for ai above. For the case of Ep/A = 200 MeV and FJq = 4.0 MeV, -Ri 
reaches the value 1.50. For Ep/A = 640 MeV where the measurement was performed and F|q 
= 4.0 MeV, we find Ri = 1.01: The gain of the total cross section due to extraband transitions 
just compensates the loss due to the reduction of the intraband cross section, resulting in a zero 
net enhancement for the total cross section compared with the harmonic limit. 

Fig. 4 shows the dependence of the cross sections and enhancement factors for the DGDR 
excitation on projectile energy. In Fig. 4(a), a2, o"™*'^^, ^'2^'^ and 0-2^*'''^ are the total, intraband, 
harmonic limit and extraband energy-integrated cross sections, respectively. (Actually, af^^'^^ is 
not a cross section itself but rather the difference of two cross sections). All of the cross sections 
increase monotonically with Ep/A. The total cross section largely stems from intraband exci- 
tation. The extraband transition contributes not more than 10% compared with the harmonic 
limit. For Ep/A < 420 MeV, o"™*^^ is larger than (72'^^. This indicates once again the damping 
enhancement which occurs at lower projectile energies where the extraband contribution is neg- 
ligible. For very low projectile energy, the enhancement factors Ri and i?™*''^ can reach several 
hundred percent. This results from the damping enhancement. When Ep/A is larger than 500 
MeV, i?2, R^^^'^^ and Rf^^'^^ are constant with values below 1.10 and 0.1, respectively. 

We turn to the dependence of the cross section for the DGDR excitation on the parameters 
FJo, F^o! FJ, F^, F^o aiid F^, using Ri and i?'/'*'''". We note that the parameters F^, F^, F^g and 
F20 appear in the expression for the cross section only in the combinations F2 + F2 and F20 + F20. 

Fig. 5(a) shows Ri and i?™*'"'^ as functions of FJg. The ratios decrease with increasing FJo, 
and are reduced by almost a factor of two when F|q is as large as several MeV. For realistic 
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Figure 3: Enhancement factors for the energy-integrated cross sections for DGDR excitation as 
described in the text as functions of FJo at different projectile energies. Parameter values are 
FJo = 0.11, F^o = 0.026, Fj = 0.30, F^ = 0.16, F^q = O-SFJo and F^ = 0.5 (MeV). 
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Figure 4: (a) Energy-integrated cross sections and (b) enhancement factors for the DGDR 
excitation versus Ep/A. Parameter values are fIq = 0.11, tIq = 0.026, FJq = 4.0, FJ = 0.30, F2 
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values of Fjo the reduction is quite small, however. An increase of FJq reduces the lifetime 
of the one-phonon state and suppresses the intraband excitation which dominates the DGDR 
excitation (Fig. 4). Both ratios depend only weakly or not at all on FJ (Fig. 5(b)). It is clear 
that the dominant intraband excitation is independent of FJ. The slow decrease of Ri with 
increasing FJ is due to the loss of extraband transition strength. However, the latter gives only 
a small contribution to -Ri. In Fig. 5(c) and Fig. 5(d) we display the dependence of Ri and of 
^mtra ^j^ "pT _j_ pi ^^^ qj^ pT^ _j_ pi^^ respectively. In both cases Ri increases with the increase 
of the variable. However, this increase is so slow that the precise choice of these parameters is 
irrelevant for our calculations. The two variable combinations on the abscissas of Figs. 5(c) and 
5(d) do not cause any loss of probability flux feeding the DGDR. Therefore, we attribute the 
slow rise of Ri in both figures to damping enhancement. This is obvious for Fig. 5(c) where the 
intraband cross section is independent of Vl + F2 and where the increase of Ri is due to the 
damping enhancement for the transition |0A;) -^ \^k). In Fig. 5(d) damping enhancement works 
also for i?™*''^. It is interesting to see that the damping enhancement survives (albeit weakly) 
even if Ep/A is as high as 640 MeV. 

We turn to a comparison of our results with experimental data [^ . In keeping with exper- 
imental results [^, we use 122% of the sum rule for all dipole matrix elements. Fig. 6 shows 
the differential cross section for the DGDR excitation as a function of excitation energy. The 
cross section has been normalized to its maximum value. The width of the DGDR is 6.1 MeV. 
Integrating the differential cross section over excitation energy from neutron threshold at 7.5 
MeV up to 40.0 MeV (this corresponds to the experimental situation), we find the value 410 
mb. This value is somewhat larger than that given in Ref. [BO], 380 ± 40 mb, where the folding 



method [^ is used. Both our DGDR width and energy-integrated cross section are consistent, 
however, with the experimental values within experimental errors. In the present case, the en- 
hancement factor Ri is 1.03, almost the same as shown in Fig. 3. In Fig. 7, we plot the ratio 
of the energy-integrated cross sections for the DGDR and the GDR as a function of Fjg. The 
ratio first increases and then decreases as this parameter increases. At the experimental value 
FJq = 4.0 MeV, the ratio is 0.0988, which again agrees with experiment (0.11 ± 0.013 ||3( 
within the experimental error. 
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Figure 5: Enhancement factors of the energy-integrated cross sections for DGDR excitation 
versus (a) FJq, (b) FJ, (c) F2 + Fg and (d) F20 + Fgg- Parameter values are Ej^j A = 640, F|q = 
4.0, F^o = O.SFio, (a) T^ = 0.026, FJ = 0.30, F^ = 0.16, F^ = 0.5; (b) T^ = 0.11, tI, = 0.026, 
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Figure 7: Ratio of energy-integrated cross sections for the DGDR and GDR excitation versus 
spreading width FIq. Parameter values are bmin = 15.5 fm, Ep/A = 640, FJq = 0.11, F 
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6 Comparison with the approach by Carlson et al. 



As mentioned in Section |T], Carlson et al. |TT], |l^ implemented the Brink-Axel hypothesis using 
a statistical model developed by Ko |]14|. The model of Ref. [ill, Q differs from the one used 
here mainly in terminology: The authors use a direct product representation of Hilbert space 
rather than the direct sum representation of our Eq. (|I]). In Refs. [|Tl|, |12] the states in Hilbert 
space are written in the form \nQ). Here n denotes the number n = 0, 1, . . . of phonons and 
a enumerates the background states. As a consequence, the mixing of the one-phonon state 
with the background states is formally described as a de-excitation of the former, i.e., as the 
transition |10) — > |0a). The Brink-Axel hypothesis is automatically implemented because dipole 
excitation affects only the first component n in the state labelled \na). As far as we can see, the 
physical content of both models is quite the same, however. The main difference between both 



approaches lies in the analytical evaluation of the average cross section. Indeed, in Refs. ||Tl| , [12 



the average over the background components labelled a is performed over the Green function 
(rather than over the intensity) and, although intuitively reasonable, does not involve controlled 
approximations. We are, in fact, not able to say whether the statistical assumptions invoked 



in Refs. [fTll , [1^] are the same as ours or not. In contradistinction to this procedure, we have 
averaged in Section ^ the intensity, and we have reduced the resulting expression by controlled 
approximations to the final form of Eq. ( |39D . We are sure that for realistic values of the decay 
widths, the statistical error of this equation and the cross section formulas derived from it, is less 
than one percent. The terminology differs also regarding the term "harmonic approximation". 
In Refs. |1TT| , [12[| , this term denotes what is here called the intraband transition. 

D are quite different in detail from the ones 



The numerical calculations of Carlson et al. ITT 



performed here. This is why we can only compare the results and not the intermediate steps. We 
do this for the reaction of ^°'^Pb on ^''^Pb at 640 MeV discussed throughout this paper. For the 
ratio of the energy-integrated cross section over the energy-integrated intraband cross section 



(we recall that this is referred to as the harmonic approximation cross section in Refs. |Tl|, p!2|), 
Carlson et al. find the value 120%. This has to be compared with our value 108%. We conclude 



that the approach of Carlson et al. [|TT| ^ constitutes a fair approximation to the exact value. 
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7 Summary 

We have studied the DGDR excitation using the Brink- Axel mechanism. The background states 
which couple to the one-phonon and two-phonon states are described in terms of the Gaussian 
Orthogonal Ensemble of random matrices. We use second-order time-dependent perturbation 
theory and calculate analytically the ensemble-averaged cross section for the DGDR excitation. 
This quantity is a function of the various decay widths and spreading widths of the one-phonon, 
two-phonon, and background states. The decay widths have been calculated from the optical 
model and the exciton model. The spreading widths are taken from experiment or used as 
parameters. We have numerically studied the dependence of the cross section for excitation of the 
DGDR on the various parameters of the theory. We have shown that for realistic values of these 
parameters, the contribution to the cross section from the Brink- Axel mechanism is significant. 
This is especially true for the reaction ^osp]^ _|_ 208p|^ ^^ projectile energy 640 MeV/nucleon. We 
have compared our results with those of Carlson et al. and have found that their work provides 
a fair approximation. For sensible values of the various spreading widths, we find that the width 
of the DGDR, the value of the energy-integrated cross section, and the ratio of this quantity over 
the energy-integrated cross section for single GDR excitation, all agree with experiment within 
experimental errors. We take this as a strong indication that the present approach accounts 
quantitatively for the DGDR excitation. It would not be difficult to use our Eq. ( |39D for the 
analysis of other data sets. Clearly, the formalism can be extended to triple-phonon excitation. 
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